Adiabatic quantum pumping through surface states in 3D topological insulators 
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We investigate adiabatic quantum pumping of Dirac fermions on the surface of a strong 3D 
topological insulator. Two different geometries are studied in detail, a normal metal - ferromagnetic 
- normal metal (NFN) junction and a ferromagnetic - normal metal - ferromagnetic (FNF) junction. 
Using a scattering matrix approach, we first calculate the tunneling conductance and then the 
adiabatically pumped current using different pumping mechanisms for both types of junctions. 
We explain the oscillatory behavior of the conductance by studying the condition for resonant 
transmission in the junctions and find that each time a new resonant mode appears in the transport 
window, the pumped current diverges. We also predict an experimentally distinguishable difference 
between the pumped current and the rectified current. 

PACS numbers: 73.20.-r, 73.40.-c, 73.23.-b, 73.63.-b 
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I. INTRODUCTION 

Recently, surface states in topological insulators have 
attracted a lot of attention in the condensed-matter com- 
munity. Both in two-dimensional (e.g., HgTe) and in 
three-dimensional (e.g., P^Sea) compounds with strong 
spin-orbit interaction the topological phase has been 
demonstrated experimentally^— . Although these com- 
pounds are insulating in the bulk (since they have an en- 
ergy gap between the conduction band and the valance 
band), their surface states support topological gapless 
excitations. In the simplest case these low-energy excita- 
tions of a strong three-dimensional topological insulator 
can be described by a single Dirac cone at the center of 
the two-dimensional Brillouin zone (r pomt)2£&~—. The 
corresponding Hamiltonian is given by^ 



Ho = hvpd ■ k — 



(1) 



Here a represents a vector whose three components are 
the three Pauli spin matrices, / represents a 2 x 2 iden- 
tity matrix in spin space, vf is the Fermi velocity, and 
\x is the chemical potential. The low-energy excitations 
of Ho [Eq. (JT}] are topologically protected against per- 
turbations^. This has prompted recent research on the 
transport properties of surface Dirac fermions. For ex- 
ample, the conductance and magnetotransport of Dirac 
fermions have been studied in normal metal - ferromag- 
net (NF), normal metal - ferromagnetic - normal metal 
(NFN) and arrays of NF junctions on the surface of a 
topological insulator— ~—, suggesting the possibility of 
an engineered magnetic switch. An anomalous magne- 
toresistance effect has been predicted in ferromagnetic- 
ferromagnetic junctions— Also, electron tunneling and 
magnetoresistance have been studied in ferromagnetic - 
normal metal - ferromagnetic (FNF) junction s 15 : 16 , for 
which it has been predicted that the conductance can 
be larger in the anti-parallel configuration of the mag- 
netizations of the two ferromagnetic regions than in the 
parallel configuration. In addition, a large research effort 
has been devoted to studying models which predict the 



existence of Majorana fermion edge states at the interface 
between superconductors and ferromagnets deposited on 
a topological insulator— : 17 : 18 . 

In this article we investigate adiabatic quantum pump- 
ing of Dirac fermions through edge states on the sur- 
face of a strong three-dimensional topological insulator. 
Quantum pumping refers to a transport mechanism in 
meso- and nanoscale devices by which a finite dc cur- 
rent is generated in the absence of an applied bias by 
periodic modulations of at least two system parameters 
(typically gate voltages or magnetic fields )^2r—. In or- 
der for electrical transport to be adiabatic, the period 
of the oscillatory driving signals has to be much longer 
than the dwell time Td we ii of the electrons in the system, 
T = 2iruj~ 1 ^> Tdweii- In the last decade, many different 
aspects of quantum pumping have been theoretically in- 
vestigated in a diverse range of nanodevices, for example 
charge and spin pumping in quantum dots22r— , the role 
of electron-electron interactions^ 6 - - — , quantum pumping 
in graphene mono- and bilayers^— as well as charge 
and spin pumping through edge states in quantum Hall 
systems^ and recently a two-dimensional topological in- 
sulator—. On the experimental side, Giazotto et aZ— 
have recently reported an experimental demonstration of 
charge pumping in an InAs nanowire embedded in a su- 
perconducting quantum interference device (SQUID). 

Our main focus is to study quantum pumping induced 
by periodic modulations of gate voltages or exchange 
fields, which are induced by a ferromagnetic strip in two 
topological insulator devices: a NFN and a FNF junc- 
tion, see Figs. Q] and [2] Using a scattering matrix ap- 
proach, we obtain analytical expressions for the angle- 
dependent pumped current in both types of junctions. 
We find that the adiabatically pumped current in a NFN 
topological insulator junction induced by periodic mod- 
ulations of gate voltages reaches maximum values at spe- 
cific energy values. In order to explain the position of 
these values, we study in detail the conductance of the 
junctions. In particular, we provide an explanation for 
resonances in the conductance that were predicted but 
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cigcnstates are given by 



FIG. 1. (Color online) Sketch of the N;FN r junction on the 
surface of a topological insulator. Pumping is induced by 
applying gate voltages (not shown) to the normal leads. In 
the middle region a thin ferromagnetic film induces ferromag- 
netism on the surface of the topological insulator by means 
of the exchange coupling^. The arrow in the middle region 
indicates the direction of the magnetization M in this region. 



not explained in detail in the previous work a 11 ' 12 . We 
show that each time a new resonant mode appears in the 
junction the conductance increases and the pumped cur- 
rent reaches a maximum value. For the FNF pump we 
predict a non-zero current by periodic modulation of the 
exchange magnetic coupling in the absence of external 
voltages. We observe and analyze basic similarities and 
differences between the two pumps studied in this paper 
and highlight an experimentally distinguishable feature 
between the pumped current and the conductance. 

The remainder of the paper is organized as follows. In 
Sec. HH we describe the NFN and FNF junctions and use 
a scattering matrix model to calculate the reflection and 
transmission coefficients of both junctions. In Sec. lIII| we 
review the conductance of the NFN junction and present 
a detailed analysis of the plateau-like steps that appear 
in the conductance. We also analyze and compare the 
conductance of the FNF junction with parallel and anti- 
parallel configuration of the magnetization. In Sec. IIV1 
we calculate the adiabatically pumped current for the 
two different pumps and derive analytical expressions as 
a function of the angle of incidence of the carriers. We 
also investigate the dependence of the pumped current 
and the conductance on the width d of the middle region. 
Finally in Sec. |V] we summarize our main results and 
propose possibilities for experimental observation of our 
predictions. 



II. NFN AND FNF JUNCTIONS 

We first describe the NFN junction, see Fig. Q] The 
junction is divided into three regions: region Ni (for 
x < 0), region N r (for x > d) and the ferromagnetic 
region F in the middle. The left and right-hand side of 
the junction represent the bare topological insulator. The 
charge carriers (surface Dirac fermions) in these regions 
are described by the Hamiltonian Hq [Eq. (JTJ] whose 
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±e ±ia 



e ±ik n x e %qy 



(2) 



where +(— ) labels the wavefunctions traveling from the 
left (right) to the right (left) of the junction. The angle 
of incidence a and the momentum k n in the a;-direction 
are given by: 



sin(a) 



hvFq 




(3) 



(4) 



Here e represents the energy measured from the Fermi en- 
ergy tp and q denotes the momentum in the y-direction. 
In the normal regions TV; and N r a dc electrical volt- 
age can be applied via metallic top gates to tune the 
chemical potential [i and thereby control the number of 
charge carriers incident on the junction. We assume gate 
voltages to be small compared to the bandgap for bulk 
states (eVi -C E g ~ 1 eV, i = Z,r), so that transport is 
well described by surface Dirac states^. In this case, the 
cigcnstates are given by 
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±e ±lai 
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D ±ik, H x iqy 



^±ik„ 



,(x-d) e iqy 



sin (ai) 



hvpq 



eVi 



' e + fi — eVi 

hv F 



q\ 



(5) 



(6) 



(7) 



(8) 



where the index i — l,r labels the normal sides of the 
junction. 

In the middle region M of the junction (0 < x < d), the 
presence of the ferromagnetic strip modifies the Hamil- 
tonian by providing an exchange field. The Hamil- 
tonian that describes the surface states is now % = 
TLq +'Hinduccd; where the induced exchange Hamiltonian 
is given b y 12 ' 14 



^induced — Hv F AI(T, 



(9) 



with the magnetization M = My. The magnitude M 
depends on the strength of the exchange coupling of the 
ferromagnetic film and can be tuned for soft ferromag- 
netic films by applying an external mag nctic fieldil The 
cigcnstates of the full Hamiltonian H are then given by: 
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(10) 
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FIG. 2. (Color online) Sketch of the F;NF r junction. Ferro- 
magnetic films are placed on top of the topological insulator 
on the left and right providing exchange fields in these re- 
gions. The arrows indicate the direction of the corresponding 
magnetizations Mi and M r , see the text for further explana- 
tion. 



with 



and 



sin(a m ) 



hv F (q + M) 
+ ' 



(11) 



Now we describe the FNF junction, see Fig. [2] Region 
Fi (x < 0) and region F r (x > d) are modeled as ferro- 
magnetic regions, respectively, with different magnetiza- 
tions Mi, M r along the y-axis and corresponding wave- 
function ipp [Eq. (fTT)|) ]. The Dirac fermions in the middle 
region N (0 < x < d) are described by the wavefunctions 
ipN [Eq.®]. When calculating transport properties of 
the FNF junction, we focus on two different alignments 
of the magnetizations of the ferromagnetic regions: the 
parallel configuration (M; || M r ), where the magnetiza- 
tions in the ferromagnetic regions point in the same di- 
rection, and the anti-parallel configuration (M; || - M r ), 
in which the magnetizations are in opposite directions. 

Using Eqns. (f2|)- (fl~2"j) we can calculate the reflection and 
transmission coefficients for a Dirac fermion with energy e 
and transverse momentum q incident from the left on the 
junction, for both the NFN and the FNF junctions. To 
this end, we consider a general F;F m F r junction, where 
the wavefunctions in each of the three regions left (7), 
middle (m) and right (r) are given by: 




-(q + Mf 



(12) 



From Eq. (fT2"j) we see that for a given energy there exists 
a critical magnetization 



M c =±2|e + ^|/(fiujr) 



(13) 



beyond which for all transverse (q) modes the wave- 
function changes from propagating to spatially decaying 
(evanescent) along the x-directio n 11 ' 12 . 



4>r 



t r l 1p+. 



(14) 



Here ipj (j = l,m,r) are the wavefunctions ([5]), ((6]) or 
(|10l) (depending on the junction considered) and ru and 
t r i denote the corresponding reflection and transmission 
coefficients. By requiring continuity of the wavefunction 
at the interfaces x = and x = d, we obtain the reflection 
and transmission coefficients: 
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(15) 
(16) 



Here ctj denotes the polar angle of the wavevector in re- 
gion j = l,m,r [Eqns. (J7|) and (|TT|) ]. When considering 
an electron incident from the right lead, one can similarly 
obtain r rr and ti r . These expressions for the reflection 
and transmission coefficients form the basis of our cal- 
culations of the conductance and the pumped current in 
Sees. [HT] and |IV] respectively. 

III. CONDUCTANCE 

The conductance G NFN of a topological insulator NFN 
junction has been studied in earlier work by Mondal et 
a h 11 ' * 12 , who predicted oscillatory behavior of G NFN as 
a function of the applied bias voltage (see also Fig. [3]). 



In this section we first briefly review their results and 
then add a quantitative explanation for the oscillations 
of the conductance. This explanation is crucial for under- 
standing the behavior of the pumped current in the next 
section. We also calculate and analyze the conductance 
in a FNF junction. 

The general expression for the conductance G across 
the junction in terms of the transmission probability 
T(a) = |t r /(a)| 2 is given by 

/■tt/2 

G=(G /2) T(a) cos a da. (17) 

J-tt/2 

Here G = ^p{eV)hv F W, p(eV) = \^ey\j (2i:(hv F ) 2 ) 
denotes the density of states, W is the sample width, 
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FIG. 3. (Color online) The conductance G NFN of the NFN 
junction [Eg. (|17[) ] as a function of e/fi for Vj = V r = 
and for different values of the effective magnetization M = 
hvpM/fi = 3 (solid blue line), 3.5 (dashed green line), 4 (dot- 
dashed red line) and 4.5 (dotted light-blue line). The effective 
junction width d = /j,d/(hvF) = 5. 



and the integration is over all the angles of incidence a. 
For Vi = V r = (ai = a r = a) the angle-dependent 
transmission probability T NFN (ct) is given b y 11 ! 12 



T 



NFN 



(a) = cos 2 (a) cos 2 (a m )/ [cos 2 (fc m d) cos 2 (a) cos 2 (a m ) 



- sin 2 (fc m d)(l — sin(a) sin(cv m ))^ 



where a m is the polar angle of the wave vector in the 
middle region as defined in Eq. HI]). This angle can be 
expressed in terms of a using the fact that the momentum 
is conserved along y-axis as: 



sin(a m ) = sin(a) + 



Mhvp 
I e + M 



(19) 



Figure [3] shows the conductance of the NFN junction 
[obtained from Eqns. (TTT]) and (|18p] as a function of the 
energy of the incoming carriers e//i for different values of 
the effective magnetization M = HvfM/h. For a given 
magnetization M, the conductance is zero for e < e c , i.e., 
G NFN (e) = 0, with the critical energy e c = hv F M/2 - fi. 
Below this energy there are no traveling modes inside the 
barrier. Our results agree with the previous results in the 
literatur o n i 12 . 

From Fig. [3] it can be observed that the conductance 
changes from plateau-like to oscillatory as e/ fi increases. 
In order to provide an explanation for this behavior we 
first analyze the plateau like regime in detail. After set- 
ting ai = a r = a in Eq. (|15p , we begin by finding the con- 
ditions when the reflection coefficient is zero, i.e., ru = 0. 
The first, trivial, condition a = a m + 27m corresponds 
to the situation of an entirely normal junction (i.e., no 
ferromagnetic region). The second and more interesting 
condition is sia{k m d) — 0. This is the case when trans- 
mission occurs via a resonant mode of the junction and 
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FIG. 4. (Color online) The transmission probability T NFN (a) 
(18) [Eq. I|18p] as a function of the angle of incidence a for differ- 
ent values of energy e/fi, (a) e/ [i — 0.7, (b) e//i = 0.9, (c) 
e/u = 1.2, (d) e/fi = 1.6, (e)e//j = 2.4, and (f) e/fi = 2.9. 
Parameters used are d = 5 and M = 3. 



can be written as (using Eqns. ((3]) and (|12)l ) 



hvf 



-d\ / 1 — I sin a 



hv F M\ _ 



(20) 



Eq. (|20|) indicates that for a given M and e there are 
certain privileged angles a c for which the barrier becomes 
transparent: 



sin(o; c ) = ± 



\ 



1 - 



M 



(21) 



with d = dfj,/(hvp) being the dimensionlcss barrier 
length. These modes are referred to as resonant modes 
in this article. 

Figure [4] shows the transmission probability [Eq. (|18l) ] 
as a function of the angle of incidence a for different 
values of energy e/fi. The dashed (red) vertical lines 
correspond to the angles satisfying Eq. (|2ip for differ- 
ent n. It can be seen that as the energy increases more 
resonant modes become available for transmission. It is 
also worth noting that for energies at which only one 
mode is present (n = 1, see Fig. QJa)) the transmission 
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FIG. 5. (Color online) The real part of the angle of incidence 
a c [Eq. (f2lj) ] versus energy e/fi for the modes n = 1, 2, 3, 4, 5. 
Parameters used are the same as in Fig. 2] 



is strongly localized at one particular angle. This prop- 
erty could be exploited to fabricate single-mode filters. 
For low energy excitations only negative angles a (i.e., 
g-momcnta anti-parallel to M) contribute to the conduc- 
tance, see Figs. |U[a)-(d). As the energy increases, the 
resonant modes move from the left to the right and also 
positive angles a (i.e., g-momenta parallel to M) begin 
to contribute, see Figs.SJe) and (f). 

Now we address the question why a mode becomes 
resonant in the barrier. Figure [5] shows the angle of inci- 
dence a c [Eq. ([21]) ] for different values of n as a function 
of energy e/fx. For a given n, the resonant mode does not 
contribute to the conductance if the energy e satisfies the 
condition e c < e < e" n , because the imaginary part of the 
momentum k m is nonzero and thus the mode is decaying 
along the x-direction. This critical energy e"™ for each 
mode is given by: 



2d 2 M 2 ' 



(22) 



When e > e"™, k m becomes real and the mode be- 
comes resonant. As we increase the energy, all the modes 
asymptotically reach their saturation angle a = n/2. 

Finally, we analyze the effect of the appearance of sub- 
sequent resonant modes on the conductance. For small 
energies, the conductance increases in plateau-like steps 
(see Fig. |3]). The first plateau corresponds to the sit- 
uation in which the first transmission mode appears at 
a = —n/2. As the energy increases, a new resonant 
mode appears and the conductance increases in a step- 
like manner. The plateaus are not sharp due to the fact 
that each new mode appearing is not sharply peaked, 
but rather has a certain distribution around a particular 
angle of incidence, see Fig. [4] Once the energy is large 
enough for there to be contributions from both positive 
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FIG. 6. (Color online) The conductance G NFN of the NFN 
junction as a function of e//j for different values of gate volt- 
ages, eV/fi = (solid blue line), 1 (dashed green line), 2 
(dashed-dot red line) and 3 (double-dotted light-blue line). 
As before, M = 3 and d = 5. 



effect of the magnetic barrier disappears and the conduc- 
tance becomes unity (G NFN = Go). 

Figure [5] shows the conductance as a function of energy 
for several values of applied bias voltages Vi = V r = V. 
As expected, the features of the conductance remain the 
same for finite V. As we increase eV, the critical energy 
e c = HvfM /2 + eV/2 — /j, for the onset of the conductance 
increases and the spacing between two consecutive reso- 
nant modes decreases. As a result the plateaus become 
narrower. 

In the remaining part of this section we study the 
conductance in a topological insulator FNF junction, as 
shown in Fig. [5] We consider both the junction with 
parallel and with anti-parallel magnetization in the fer- 
romagnetic regions. In the parallel configuration, us- 
ing ai = a r = a = sin - {%vp{q + M)/|e + fi\) and 
a m = sin -1 (Hvpq/\e + fi\) in Eq. (I16|) . we find that the 
conductance is similar to the conductance of a NFN junc- 
tion, as displayed in Fig. [3] However, in the FNF junction 
the first resonant mode becomes resonant for positive a 
(i.e., transverse q- momentum parallel to M) and as the 
energy increases, the resonances move towards negative 
values of the angle a. This, however, does not affect the 
total conductance, as we sum over all possible angles of 
incidence, and the same analysis as for the NFN junction 
presented above can be applied to understand the FNF 
junction with parallel magnetization. 

In the case of anti-parallel alignment of the magne- 
tization in the two ferromagnetic regions, we substitute 
sin(az) = HvF{q+M)/\e+[i\, sin(a r ) = hvp{q— M)/\e+fj,\ 
and sin(a m ) = hvpq/\e + jj,\ in Eq. (|16p . The conduc- 
tance of this junction was studied previously in Rcfs i 15 ' 16 
and the transmission probability T FNF ' AP (ai,a r ) = 
\t r i(ai,a r )\ 2 is given by: 
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rpFNF.AP/ s COS 2 (o;,)cOS 2 (a m ) 

1 ' (ai,a r j = o. (23) 

cos 2 (fc m d) cos 2 (^i±^)cos 2 (a m ) + sin 2 (fc m d) [cos (^f^) - sin(2i±^) sin(a m )] 




FIG. 7. (Color online) The conductance G* FNF ' AP of the FNF 
junction in the anti-parallel configuration as a function of e//j,, 
for M = 3 (blue solid line), 3.5 (green dashed line) and 4 (red 
dotted line). The parameter d = 5. 



The total conductance is obtained by multiplying 
T FNF,AP (a;, a r ) with cos(av)/ cos(a/) and then integrat- 
ing over the allowed angles of incidence^, i.e., from a ci = 
sin~ 1 (2hv F M/(\e + fi\)-l) to a C2 = sin^ 1 (2hv F M/(\e + 
fi\) + 1). Thus we can write 



g fnf,ap = Go/2 / G tNt ' At, ( ai , ar) cos{ai)dai 



tt/2 



FNF,AP / 



where 



,oFNF,AP/ „ \ CQs(q r ) FNF ap 

G (ai,a r ) = ; — rJ 



(24) 



(a h a r ). (25) 



cos(a/) 

Fig. [7] shows the conductance of the FNF junction in 
the anti-parallel configuration. From the horizontal axis 
we see that the critical energy e c for the onset of the 
conductance is larger than in the corresponding parallel 
configuration. Moreover, as the energy e/fi increases the 
conductance exhibits no plateau behavior: it increases in 
an oscillatory fashion. This oscillatory behavior can be 
understood from Fig. [5J which shows G FNF ' AP (a) for four 
different values of e/fi. Note that all the angles (a/, a m 
and a r ) can be expressed in terms of one angle, which we 
choose to be ai = a. As the energy increases, the area 
under the curve oscillates resulting in oscillations in the 
conductance. 

Summarizing, we have obtained a quantitative expla- 
nation for the behavior of the conductance in topological 
insulator NFN and FNF-j unctions in terms of the num- 
ber of resonant modes in the junction. This explanation 
forms the basis for understanding the behavior of the 
pumped current in the next section. 
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The angle-dependent total transmission 
T" " r ,At " (a;,Qr) for the anti-parallel configuration of 
the FNF junction as a function of the angle of incidence ai 
for (a) e/fx = 2.47, (b) e/fi = 2.57, (c) e/fj, = 3.10 and (d) 
e/fi = 3.40. Parameters used are M = 3 and d = 5. 



IV. ADIABATICALLY PUMPED CURRENT 

In this section we investigate adiabatically pumped 
currents through NFN and FNF junctions in a topo- 
logical insulator—. In general, a pumped current is 
generated by slow variation of two system parameters 
X\ and X 2 in the absence of a bias voltagoJ^^. For 
periodic modulations Xi(t) = X\^ + 6X1 cos(ujt) and 
X 2 (t) = X 2i o + 5X 2 cos(uit + 0), the pumped current I p 
into the left lead of the junction can be expressed in terms 
of the area A enclosed by the contour that is traced out 
in (X\, A^-parametcr space during one pumping cycle 
as 



,20. 



L p = ^2 I dX ^ dX * Y. n (^i^2) (26a) 
SX 1 SX 2 sin<£ n Pfi>^2), (26b) 



we 
2^ 



with 



U{X 1 ,X 2 )=lm 



dr* u dr u 
dX l dX 2 ' dX 1 dX 2 



dt* r dti r 



(27) 



Here ru and U r represent the reflection and transmission 
coefficients into the left lead and the index m sums over 
all modes (a similar expression can be obtained for the 
pumped current into the right lead). Eq. (|26b[) is valid 
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in the bilinear response regime where SXi <C -X"i,o and 
6X 2 <C X 2 .o and the integral in Eq. ()26a|) becomes inde- 
pendent of the pumping contour. 

First we analyze the NFN pump, where the pumped 
current is generated by adiabatic variation of gate volt- 
ages Vi and V r which change the chemical potential in 
the normal leads on the left and right of the junction, re- 



spectively (see Fig.QJ. Calculating the derivatives of the 
reflection and transmission coefficients r« [Eq. (|T5|) ] and 
t r i with respect to a/ and a r , substituting into Eq. (|27j) 
and using dotj/(edVj) = tan(aj)/|e + fj, — eVj\ (j = Z,r), 
the pumped current for Vi = V2 = V and for a specific 
angle of incidence a is given by: 



jNFN/ x jNFN cos 3 (a m ) sin 2 (q) cos(a) sin(2fc m d) 

(1 + e/jj, — eV j /j,) 2 (cos 2 (a) cos 2 (a m ) cos 2 (k m d) + sin 2 (fc m d)(l — sin(a) sin(a m )) 2 ) 2 
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FIG. 9. The pumped current 7^ FN [Eq. 
junction as a function of e/ fj, for V = 0, d 
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p0|)] in the NFN 

-- 5 and M = 3. 



Here I$ bN = we/(87r) smty) (eSVi / fj.)(e6V 2 / fi) and 
sin(a m ) is given by Eq. (flUf. In the limit M — >• (i.e., 
a m — > a) in an entirely normal junction, we obtain from 
Eq. (|2"g)l the angle- dependent pumped current as: 



tNFNi 

1 p \M=0 



t-NFN 
J 



tt/2 



tt/2 



i(2A; m d) 



sm a 



(1 + e//x — eV/fj,) cos 3 a 



da. 



(29) 

On the other hand, the transmission T NFN (a) [Eq. (f!8| ] 
in this limit is given by T NFN | fl7 ^ -> 1, independent 
of the angle of incidence a. We notice that even if the 
probability for transmission is one, it is possible to pump 
a current in the adiabatic driving regime. The total 
pumped current I^ FN is then obtained by integrating 



over a: 



rNFN 



tt/2 



-tt/2 



Ip FN (a) cos a da. 



(30) 



In general, this integral cannot be evaluated analytically 



and we have obtained our results numerically. Figure [9] 
shows the total pumped current /^ FN (in units of Iq FN ) 

at zero bias V = V\ = V<i = for M = 3. Comparing 
Figs. [3] and [9] we see that there is a correlation between 
the pumped current and the conductance for the NFN 
junction: for low energies, the pumped current I^ FN is 
zero as no traveling modes are allowed in the junction. 
As we increase the energy, each time a resonant mode 
appears [see Eq. ([20)) ]. the pumped current diverges and 
changes sign. For energies where both positive and neg- 
ative angles of incidence contribute to the conductance, 
the pumped current remains finite but keeps changing 
its sign. From Figs. [3] and [9] it can also be seen that 
the pumped current vanishes for energies at which sub- 
sequent resonant modes become fully transmitting. In 
order to gain further insight we plot the analogue of Fig. @] 
for the pumped current. Figure 1101 shows the pumped 
current [Eq. [3D] as a function of the angle of incidence a 
for different values of e//Lt. The chosen values of e//i are 
same as in Fig. 0] We see that the features in Fig.fTUlhavc 
a direct correlation with the features in Fig. [4] whenever 
there is a sharp peak in the transmission the pumped 
current diverges and changes sign. The key feature that 
distinguishes between the pumped current and the con- 
ductance is that the pumped current changes sign at par- 
ticular values of the energy, while the conductance does 
not. 



Now wc analyze the pumped current in the FNF junc- 
tion with parallel orientation of the magnetizations. In 
this system, the driving parameters are the magneti- 
zations Mi and M r in the left and right contacts, re- 
spectively, see Fig. [5[ After calculating the derivatives 
of the reflection and transmission coefficients, obtaining 
the imaginary part of Eq. ([27|) . and using daj/dMj = 
hvp/{\e + fx\ cos(aj)) (j = l,r), the pumped current 
(a) for Mi = M 2 = M is: 



7FNF 

V 



rFNF / \ rFNF cos (a m ) cos(a) sin(2/c m d) 

1 (a) = — 1q 5 . (31) 

(1 + e / /i) 2 (cos 2 (a) cos 2 (a m ) cos 2 (fc m d) + sin (fc m d)(l — sin(a) sin(a m )) 2 ) 2 
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■7/2 1 -jr/4 



(a) 




jr/4 



jr/4 



x/4 



?r/2 



(=) 



7T/2 



(e) 



t/2 




E -0.4 




jr/4 )r/2 



(</) 



ir/4 



(M locations where the conductance changes sharply. But 
there is an important difference between both pumped 
currents. The pumped current in an NFN-junction at 
normal incidence vanishes, Ip™(a = 0) = 0, while 
/ FNF (a = 0) 0. This difference arises because the 
two pumps are driven by two different parameters (volt- 
ages in the NFN pump and magnetizations in the FNF 
pump). 

Finally, we briefly analyze the behavior of the pumped 
ir/2 current as a function of the width d of the middle region. 
For energies below e c , the pumped current of the NFN 
junction decays to zero as the width d increases (there 
are no resonant modes in the system). For energies larger 
than e > e c , the pumped current i^ FN oscillates as a 

function of width d. Fig. [TT] shows the conductance and 
the pumped current as a function of d for e//j, = 2.5. The 
peaks in the conductance correspond to the resonance 
ir/4 tt/2 condition Eq. (|2T)|) . The pumped current i^ FN changes 

its sign at exactly the same values of d where the con- 
ductance has a maximum. This analysis holds as well for 
the FNF junction. 



(/) 



FIG. 10. The pumped current Ip FN (a) [Eq. |ggj] as a func- 
tion of the angle of incidence a for different values of energy 
e/n, (a) e/n = 0.7, (b) e/fi = 0.9, (c) e//i = 1.2, (d) e//i = 1.6, 
(e) e/fj, = 2.4, and (f) e/fi = 2.9. Parameters used are d = 5 
and M = 3. 





FIG. 11. (a) The conductance of the NFN junction as a 
function of d. (b) The pumped current for the NFN junction 
as a function of d. Parameters used are e/fj, = 2.5, M — 3 
and V = 0. 



Here 7 FNF = ue/(8n)sm(<t>)(hv F 6Mi/fi)(hv F 6M r /pi) 
and sin(a m ) = sin(a) — hv F M/(\e + fi\). 

The behavior of the current J FNF is similar to that 
of the pumped current i^ FN in a NFN-junction (shown 
in Fig. [9]). This can also be seen by comparing the de- 
nominators in Eqns. (|28|) and ([31]) . Again we observe 
that the pumped current diverges at exactly the same 



V. SUMMARY AND DISCUSSION 

To summarize, we have analyzed quantum transport 
by Dirac fermion surface states in NFN and the FNF 
junctions in a 3D topological insulator. We have shown 
that for low energies the appearance of a new resonant 
mode results in a plateau-like increment of the conduc- 
tance and a diverging pumped current in these junctions 
which also changes sign. This is our key result, and rep- 
resents an experimentally distinguishable signature be- 
tween conductance and the pumped current. We high- 
lighted an interesting difference between the two differ- 
ent pumping mechanisms for the NFN and FNF junc- 
tions, observing different behaviors for normal incidence 
(q = 0). Experimentally, the NFN pump could be re- 
alized using current technology The FNF pump will be 
more difficult to realize since it requires oscillating mag- 
netizations. A possible way to realize a FNF pump could 
be by moving the two ferromagnetic layers coherently 
using a nanomcchanical oscillator—. Experimental veri- 
fication of our predictions will provide further insight into 
quantum transport through these junctions. 
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